Mg/Ti multilayers: structural, optical and hydrogen absorption properties 
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Abstract 

Mg-Ti alloys have uncommon optical and hydrogen absorbing properties, originating from a "spinodal-like" microstructure 
with a small degree of chemical short-range order in the atoms distribution. In the present study we artificially engineer short- 
range order by depositing Pd-capped Mg/Ti multilayers with different periodicities and characterize them both structurally 
and optically. Notwithstanding the large lattice parameter mismatch between Mg and Ti, the as-deposited metallic multilayers 
show good structural coherence. Upon exposure to H2 gas a two-step hydrogenation process occurs, with the Ti layers forming 
the hydride before Mg. From in-situ measurements of the bilayer thickness A at different hydrogen pressures, we observe 
large out-of-plane expansions of the Mg and Ti layers upon hydrogenation, indicating strong plastic deformations in the films 
and a consequent shortening of the coherence length. Upon unloading at room temperature in air, hydrogen atoms remain 
trapped in the Ti layers due to kinetic constraints. Such loading/unloading sequence can be explained in terms of the different 
thermodynamic properties of hydrogen in Mg and Ti, as shown by diffusion calculations on a model multilayered systems. 
Absorption isotherms measured by hydrogenography can be interpreted as a result of the elastic clamping arising from strongly 
bonded Mg/Pd and broken Mg/Ti interfaces. 

PACS numbers: 68.65.Ac, 88.30.rd 

Keywords: magnesium, titanium, multilayers, hydrogen 
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I. INTRODUCTION 

Magnesium and titanium are immiscible. Metastable 
Mg-Ti alloys have nevertheless been successfully pre- 
pared in thin films by high energy processes, such as 
electron beam deposition 1 - and magnetron sputtering^i 3 - 
These films have gravimetric hydrogen storage capaci- 
ties up to 6.5 wt%i and fast and reversible kinetics of 
hydrogen absorption and dcsorptioni^ The structural re- 
versibility of Mg-Ti alloys is particularly surprising when 
considering the strong segregation occurring in many Mg- 
based binary systems4£ Furthermore, when exposed to 
H2, Pd-capped Mg-Ti thin films switch from a reflect- 
ing metallic state to a black, light-absorbing, hydro- 
genated stated This reversible optical black state can 
be applied in hydrogen sensors^ and smart absorbers 
for solar collectors ^ As also suggested by first princi- 
ple calculations^ the black appearance of the hydride is 
due to the formation of a face-centered-cubic phase, in 
which Mg and Ti atoms are distributed among the lat- 
tice sites with a certain degree of chemical short-range 
order^ 1 ^ leading to the coexistence of Mg-rich and Ti- 
rich regions with structurally coherent boundaries. Such 
"spinodal-like" microstructure is not uncommon in im- 
miscible binary alloy a 11 ! 12 and it is a key ingredient in 



understanding the exceptional reversibility of Mg-Ti thin 
films i 10 ' 13 In order to achieve a deeper understanding of 
the role of chemical segregation on the structural, optical 
and hydrogen absorbing properties of these systems, we 
engineered one dimensional short-range order by deposit- 
ing several Mg/Ti multilayers with different periodicities. 
By means of optical and structural studies we are able 
to reconstruct the hydrogen loading sequence, measure 
the out-of-plane expansion of the individual layers and 
detect the breaking of structural coherence occurring at 
the Mg/Ti interfaces upon formation of TiH2. Such re- 
moval of coherence is responsible for the "scissor" effect 
observed in Ti-sandwiched Mg thin filmsi when a thin 
Mg film is sandwiched between Ti layers it absorbs hy- 
drogen at pressures close to bulk Mg, effectively behav- 
ing as quasifrec. This is due to the fact that Mg and Ti, 
thanks to their positive enthalpy of mixing, form inter- 
faces with poor adhesion, which become even more dis- 
connected when Ti expands upon hydrogen absorption. 
On the contrary, in Mg films capped with Pd, alloying 
occurs at the Mg/Pd interface and Mg feels an clastic 
constraint due to the presence of the cap layer that leads 
to higher equilibrium pressures of hydrogen absorption*^ 
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II. EXPERIMENTAL DETAILS 

Mg/Ti multilayers covered with Pd are deposited in a 
ultra-high vacuum (UHV) compatible system (base pres- 
sure = 1CP 6 Pa) equipped with a computer controlled 
shutter system, by DC and RF magnetron sputtering of 
Mg (99.95%), Ti (99.999%) and Pd (99.98%) targets in 
0.3 Pa of Ar, on substrates kept at room temperature. 
The substrates used are 10x10x1 mm polished single- 
crystal Si(100) for XRD, XRR and HRTEM, silicon ni- 
tride membranes for in-plane TEM, 20x10x1 mm quartz 
for optical spectroscopy and 10x10x0.5 mm float glass for 
hydrogenography measurements. In order to obtain ho- 
mogenous films the substrates are continuously rotated 
during sputtering. The films are covered with Pd to pre- 
vent oxidation and promote hydrogen dissociation and 
absorption. The Pd thickness varies between 1 and 10 
nm depending on the experimental technique used to an- 
alyze the samples. Typical deposition rates are 0.22 nm/s 
for Mg at 150 W (RF), 0.08 nm/s for Ti at 200 W (DC) 
and 0.11 nm/s for Pd at 50 W (DC). In each Mg/Ti sam- 
ple the Mg layers are twice as thick as the Ti ones, giv- 
ing, after correcting for the molar volumes (Vyig = 13.97 
cm 3 /mol, Vti = 10.64 cm 3 /mol), an overall composition 
of Mgo.6oTio.4o- The total multilayer thickness is 60 nm 
with 6 different bi-layer thicknesses A, ranging from 1.5 
to 60 nm: Nx[Ti(§ nm)Mg(f nm)], with N = 40, 20, 
10, 5, 2 and 1. These samples characteristics allow us to 
have fast and comparable kinetics of hydrogen absorp- 
tion and desorption in all the multilayers. A sketch of 
the samples geometry is shown in Fig. [TJ 




FIG. 1: (color online) Sample geometry. The total multilayers 
thickness is 60 nm and, in every sample, Mg layers are twice 
as thick as Ti ones: Nx[Ti(f2 nm)Mg(^f nm)], with N = 40, 
20, 10, 5, 2 and 1. 



XRD patterns are measured in a 9 — 26 configura- 
tion, with a Bruker D8 Discover diffractometer equipped 



with a two-dimensional detector for real-time data col- 
lection over a large area with high sensitivity and low 
background. A beryllium dome allows in-situ diffraction 
measurements during hydrogenation/deydrogenation of 
the films in hydrogen pressures up to 10 5 Pa and tem- 
peratures between room temperature and 473 K. 

XRR measurements are performed on the lOx sam- 
ple in a specially designed UHV chamber mounted on a 
Bruker Discover D8 X-ray diffractometer equipped with 
a parallel X-ray beam (CuKai A = 0.15406 nm), that 
allows in-situ hydrogen loading in a wide range of tem- 
peratures and pressures^ The sample is measured at 
333 K both in the as-deposited metallic state and during 
hydrogen uptake at different H2 pressures. The tempera- 
ture is high enough to promote fast kinetics of hydrogen 
absorption but low enough to avoid severe alloying at the 
Mg/Pd interface. In-plane resistance measurements are 
used to determine whether the dissolved hydrogen is in 
equilibrium with the surrounding H2 atmosphere. 

HRTEM was performed on a JEOL 2100F FEG-TEM, 
equipped with Digital STEM, Gatan Orius imaging sys- 
tem, and Gatan DigiScan. Cross section samples were 
prepared by ion beam milling using an FEI Quanta 
200 ED FIB-SEM and transferred to the TEM using a 
Gatan HHST 4004 environmental cell and vacuum trans- 
fer holder to minimise oxidation. 

Hydrogen loading isotherms are measured at 333 K 
by means of Hydrogenography,— an optical technique 
that allows to detect the amount of light transmitted 
through a thin film, while slowly increasing the hydro- 
gen pressure at constant temperature. The Pressure- 
Optical transmission-Isotherms (PTIs) obtained by hy- 
drogenography can be directly related to the stan- 
dard Pressure-Composition-Isotherms (PCIs) measured 
for metal hydrides.— Details of the hydrogenography ex- 
perimental setup can be found in Gremaud et al..— 

Optical spectra are measured with a Perkin Elmer 
Lambda 900 diffraction grating spectrometer (0.5 < 
hut < 6.5 eV). Reflection and transmission spectra of 
the as-deposited and hydrogenated films are measured 
through the transparent quartz substrate at 333 K in 
vacuum and in 10 5 Pa H2, respectively. 



III. RESULTS AND DISCUSSION 

A. Structural characterization: XRD, XRR and 
TEM 

1. As-deposited multilayers 

The uncorrected diffraction patterns measured for the 
as-deposited samples in vacuum at room temperature are 
shown in Fig. [21 Both Mg and Ti have hexagonal-closed- 
packed structures in their elemental form and, as will 
be shown below, the Mg/Ti multilayers grow with the 
(001) plane parallel to the substrate surface. For the lx 
sample the Ti and Mg layers are thick enough to give 
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FIG. 2: (color online) Measured (dots) and simulated (lines) 
XRD patterns of Nx[Ti(22 nm)Mg(|2 nm )] multilayers with 
N = 40, 20, 10, 5, 2 and 1, covered with 5 nm of Pd. The 
simulations are based on idealized step model, not capturing 
the influence of defects and other imperfections. 



the reflections of pure elements, while already for the 
5x sample superlattice peaks (satellite peaks) begin to 
appear, due to the repetitions of the bi-layer thickness 
A. The broad peak appearing at ~40° for the 40x sam- 
ple is the (111) reflection from the face-centered cubic 
Pd. The absence of any Pd reflection in the patterns 
of the other multilayers is not surprising when consid- 
ering that Pd is deposited on top of the uppermost Mg 
layer. Pd deposited on Mgj,Tii_ y thin films, in fact, does 
not produce any diffraction signal for y > 0.9, although 
its presence is confirmed by Rutherford Backscattering 
Spectrometry (RBS) This peak disappearance is most 
likely due to the increase in lattice mismatch between 
Pd and Mgj,Tii_j, with increasing magnesium content, 



suggesting that for the 40x sample the Ti and Mg layers 
are thin enough to undergo lattice deformations similar 
to the ones expected in a Mgo.6Tio.4 alloy. In Fig. [2] 
we also show the simulated patterns obtained with an 
ideal step model, which assumes perfect superlattices 
with a square-wave composition modulation along the 
growth direction, and coherent Mg/Ti interfaces . 19 ! 20 In 
the model the diffracted intensity / is given by: 
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where 7n is a term due to the N bilayer repetitions in the 
multilayers, Tjyig and 7xi ar e the intensities of the con- 
situent materials and /MgTi is a mixed term arising from 
the structural coherence.— In incoherent multilayers the 
mixed term disappears leading to a different distribution 
of intensities in the satellite peaks. The scattering pow- 
ers of the elements are approximated with their atomic 
numbers and the patterns are filtered with a Gaussian 
distribution, with a full width at half maximum (fwhm) 
of 0.47°, to take into account the instrumental broad- 
ening and the possible deviations from a perfect geom- 
etry. These deviations include random variations in the 
number of atomic layers, interface roughness or interdif- 
fusion, variations in the lattice spacings due to in-plane 
clastic coherency strain, distribution of sizes in in-plane 
grains and random orientations of the growth direction 
of each grain. Several models have been developed to 
include corrections to the ideal step model and take into 
account deviations from a perfect geometry. 19 i 21 In our 
XRD measurements, however, the experimental broad- 
ening is very high, of the order of 0.3°, due to the use of 
a two-dimensional detector that has a finite grid and re- 
quires a rather large beam spot on the sample in order to 
produce a high signal-to-noisc ratio. Such a large exper- 
imental broadening hinders any attempt to obtain more 
detailed informations from the XRD patterns. Neverthe- 
less, the simulations shown in Fig. [2] qualitatively repro- 
duce all the features observed experimentally, showing 
that upon deposition Mg and Ti form well defined layers 
with partially coherent interfaces. The structural coher- 
ence in Mg/Ti multilayers is rather surprising given the 
8.7% lattice mismatch between the in-plane cell parame- 
ter of Mg (a hcp = 0.3209 nm) and Ti (a hcp = 0.2951 nm). 
The rocking curve measured over the fundamental peak 
(s = 0) of the lOx sample has a fwhm of 5.7°, indicating 
a moderately textured multilayer. 

For a multilayer the Bragg law can be written as>^ 
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where s is the index of the satellite peak with respect to 
the fundamental peak (s = 0), as shown for samples 20x, 
lOx and 5x in Fig. [Jl ^Mg/Ti is the average interplanar 
distance of a Mg/Ti bilayer and A is the bilayer thick- 
ness. In Fig. [3] we plot the ratio sin 9 s / AcuKai versus 
s, for the satellite peaks of the 20x, lOx and 5x samples. 
Using Eq. [21 from the fitted slopes, S, and intercepts, 
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FIG. 3: (color online) Measured linear dependence of the ratio 
sin# s /AcuKai on the satellite index s for the 20x (circles), lOx 
(squares) and 5x (triangles) samples as indicated in Fig. [2] 
The dashed lines are linear fits to the data with slopes S and 
intercepts I. 



I, we obtain the bilayer thicknesses A and the average 
interplanar distance dMg/Tb see Table HI The excellent 



TABLE I: Comparison between the experimental and nominal 
values for the bilayer thickness A and the average interplanar 



distance d 
5x. 



Mg/Ti 



obtained from eq. [5] for samples 20x, lOx and 



Sample 


A(nom) 


A(exp) 


4i E /Ti(nom) 


<4i g /Ti(exp) 




nm 


nm 


nm 


nm 


20x 


3.0 


3.0 


0.2518 


0.2504 


lOx 


6.0 


6.0 


0.2518 


0.2505 


5x 


12.0 


11.9 


0.2518 


0.2508 



agreement between experiment and model indicates that 
the textured regions of our Mg/Ti multilayers correspond 
closely to their ideal geometries. 

A rough estimate of the average coherence length £, de- 
fined as the distance over which the atomic positions are 
quantitatively correlated^ can be obtained with Scher- 
rer's formula: £ = (K ■ Ac U Kc«i)/(/3 • cos 9), where K is 
the shape factor (K ranging between 0.9 and 1, depend- 
ing on the grains shape) and (3 is the fwhm expressed in 
radians. For the fundamental peak (s = 0) of the lOx 
sample we have /3 w 0.5° = 8.7- 10" 3 rad and 6 = 17.92°, 
corresponding to a coherence length £ « 17 nm, assum- 
ing spherical grains. Given the very high experimental 
broadening in our XRD measurements, this value has to 
be taken only as a lower limit. Considering that for the 
lOx sample the bilayer thickness is A = 6 nm we can con- 
clude that, in the textured regions of the multilayer, the 
crystal registry is maintained at least for few Mg/Ti bi- 
layer repetitions and that the Mg/Ti interfaces are there- 
fore partially coherent, as already suggested by the com- 



parison between experimental and simulated XRD pat- 
terns (Fig. (2). 

The lOx sample, 10x[Ti(2 nm)Mg(4 nm)], has been 
explored in more detail by means of in-situ X-Ray Re- 
flectivity at different hydrogen pressures. The XRR mea- 
surement of the as-deposited state of sample lOx, covered 
with 10 nm of Pd, is shown in Fig. the measurement 
is conducted in vacuum (base pressure = 10~ 7 Pa) at 
333 K and the sample is deposited on a Si substrate with 
(100) orientation. From the satellite peaks positions in 
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FIG. 4: (color online) (a) X-Ray Reflectivity pattern of the 
lOx sample, 10x[Ti(2 nm)Mg(4 nm)]Pd(10 nm), taken at 333 
K in vacuum (base pressure = 10 -7 Pa). The red solid line is 
a simulation of the multilayer using GenX.— The inset shows 
the XRD pattern measured at high angles on a lOx sample 
deposited in the same run. (b) Real part of the Scattering 
Length Density (SLD) profile corresponding to the simulation 
in Fig. UK- 



the XRD pattern measured at high angles on a sample 
deposited in the same run (inset in Fig. [4^,) we obtain 
A = 5.63 nm and ^Mg/Ti = 0.2500 nm. The XRR curve 
is simulated with GenX-22 using the model in Fig. [JJ to 
obtain the individual thicknesses of the constituent layers 
and the roughnesses at the interfaces. The software uses 
a dynamic optical model that incorporates effects arising 
from refraction, x-ray absorption, multiple scattering, in- 
strumental resolution and instrumental geometry. In the 
fit all the layer thicknesses and the interface roughnesses 
are varied. A layer of 1.7 nm of Si02 is added to the 
simulation on top of the Si (100) substrate to take into 
account the substrate surface oxidation. The results for 
the as-deposited state of sample lOx are summarized in 
the third column of Table [Til The agreement between 
the A values obtained from the satellite peaks position 
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TABLE II: Nominal and measured structural parameters of 
the lOx sample, 10x[Ti(2 nm)Mg(4 nm)]Pd(10 nm), as ob- 
tained from fitting the XRR patterns measured in vacuum 
(as-deposited) and in 10 3 Pa of H2 (fully hydrogenated): bi- 
layer thickness A, layers thickness d and roughness a. 
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XRR 


XRR 






(vacuum) 


(Ph 2 = 10 raj 




nm 


nm 


nm 


A 
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5.6 


6.9 


rfsi0 2 




1.7 


2.0 


d-Ti 
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2.0 


2.3 


^Mg 
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3.6 


4.6 


dpd 


10 


8.7 


8.7 


°"Ti/Mg 





0.74 


1.0 


°"Mg/Ti 





0.55 


1.2 


°"Mg/Pd 





1.2 


1.3 


^Pd/vacuum 





0.78 


1.1 



and from the simulation of the XRR pattern is excellent. 
The layers are flat with interface roughnesses of the or- 
der of one unit cell (assuming hep Mg with c = 0.521 nm 
and hep Ti with c = 0.469 nm). Fig. HJd shows the real 
part of the Scattering Length Density (SLD) profile cor- 
responding to the simulation in Fig. fjji. The SLD, which 
is given by the mass density profile times the scattering 
lengths of the elements, gives an idea of the deviations 
from a perfect square-wave model that we have to intro- 
duce in the simulation, in order to accurately reproduce 
the experimental measurement. 

Another method to determine the bilayer thickness A, 
in the limit of a kinematic approximation, is by looking 
at the position of the reflectivity peaks at low angles. We 
have>2£ 

sin 2 0„ = (^£^1-^^ + 26 (3) 

where n is the order of the reflectivity peaks, as shown in 
Fig. 0k, and S is the deviation from unity of the real part 
of the average refractive index and, in first approxima- 
tion, can be neglected i 16 ' 21 From a linear fit of the plot 
sin 2 S n / X^ uKai versus n 2 we calculate A = 5.64 nm, in 
agreement with the values obtained both from the XRR 
simulation and from the satellite peak positions in the 
XRD pattern. 

In Fig. [3Ji we show the cross-section bright-field Scan- 
ning Transmission Electron Microscopy (STEM) image 
of a 20x[Ti(2 nm)Mg(4 nm)] multilayer, deposited on a 
Si(100) substrate and covered with 10 nm of Pd. The 
sample is well layered with slightly wavy Mg/Ti inter- 
faces. Due to cumulative roughness, only 16 out of 20 
bilayer repetitions are visible. Figurc[5j3 shows the inten- 
sity profile of the area delimited by the white dashed line 
in Fig. [5^. The maxima in the profile correspond to the 
bright Mg layers, except for the first peak on the bottom, 
which is attributed to the SiC>2 film covering the Si sub- 
strate. The effect of cumulative roughness is clearly visi- 
ble in the intensity profile, where the peak-to- valley ratio 
decreases with increasing film thickness. Note that the 
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FIG. 5: (a) Cross-section bright-field STEM image of a 
20x[Ti(2 nm)Mg(4 nm)] multilayer deposited on a Si(100) 
substrate and covered with 10 nm of Pd. (b) Intensity profile 
of the area delimited by a white dashed line. The maxima 
correspond to the Mg bright layers, except for the first (bot- 
tom) peak, which is attributed to Si02 at the substrate/film 
interface. 



sample measured in TEM is twice as thick as all the other 
investigated samples, in which the cumulative roughness 
effect is therefore going to be much smaller. Figure [5^ 
shows a cross-section image of the same sample, in which 
it is clearly visible how the film is partially crystalline, 
with grains extending for several Mg/Ti repetitions, as 
already suggested by X-Ray Diffraction results. Figures 
HJd and[6b show a High Resolution TEM (HRTEM) image 
of the 20x[Ti(2 nm)Mg(4 nm)] multilayer and its Fourier 
transform, respectively. As expected, the multilayer has 
hexagonal closed-packed structure and grows along to the 
[001] direction. Although a slight decrease in crystallinity 
is observed with increasing thickness (not shown here), 
the crystal structure and orientation shown in Figure [SJd 
are visible across the multilayer from the substrate up to 
the Pd cap. 



2. Hydrogen loading 

When exposed to H2 gas at room temperature, the 
Pd-cappcd Mg/Ti multilayers hydrogenate in two con- 
secutive steps: (i) at lower H2 pressures only the Ti layers 
form a hydride while Mg remains in its metallic state, (ii) 
at higher H2 pressures also Mg absorbs hydrogen form- 
ing MgH2. Such loading sequence is due to the lower 
(more negative) enthalpy of hydride formation of TiH2 
(-65 kJ(mol H)" 1 )^ with respect to MgH 2 (-37.2 kJ(mol 
H) -1 ) 2 ^ and can be detected from XRD, XRR and opti- 
cal spectroscopy measurements. 

Figure [7] exhibits the XRD patterns measured during 
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FIG. 6: (a) Cross-section TEM image of the 20x sample shown 
in Fig. [5] The dashed area highlights a region in which struc- 
tural coherence is maintained across multiple Mg/Ti inter- 
faces, (b) HRTEM image of sample 20x, taken approximately 
in the middle of the multilayer. The dashed arrow, parallel to 
the [001] direction of the hep structure, indicates the growth 
direction of the multilayer, (c) Fourier transform of Fig. [6jz>. 



loading and unloading of sample lOx. Upon loading the 
sample is exposed to 1 bar H2 pressure at room tem- 
perature. The loading sequence shows an intermediate 
step in which only one peak is visible. This intermediate 
peak is similar to what is expected for a perfect hep- (002) 
Mg / fcc-(lll) TiH 2 multilayer, (dashed lines in Fig. EJ), 
although it is shifted to slightly lower angles, suggesting 
that Mg layers might also be partially hydrogenated. The 
final hydrogenated state shows poor crystallinity with a 
broad reflection at 29 rj 34° . When unloading the multi- 
layer in air at room temperature the intermediate peak is 
recovered, suggesting that, due to kinetic limitations, hy- 
drogen does not desorb from the Ti layers in the sample. 
In order to recover the initial metallic state the sample 
has to be heated in air at 433 K. Although the satel- 
lite peaks in the final de-hydrogenated state are broader 
than in the as-deposited initial sample, it is remarkable 
that crystallinity appears again, after the large changes 
in volume upon hydrogenation and dchydrogenation of 
the Mg and Ti layers. 

An intermediate state, similar to the one observed 
upon desorption in air at room temperature for the lOx 
sample, corresponding to a hcp-(002) Mg / fcc-(lll) 
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FIG. 7: XRD patterns measured during loading (left) and 
unloading (right) of the lOx sample, 10x[Ti(2 nm)Mg(4 
nm)]Pd(10 nm). The intensities are normalized by the mea- 
suring time. The dashed lines are simulations of a perfect 
hcp-(002) Mg / fcc-(lll) TiH 2 lOx multilayer. 



TiH2 multilayer, is obtained for the desorbed state of 
all the multilayers explored in the present work, Fig. [5] 
The lines are simulations based on the ideal step model, 
in which we assume that the samples have an ideal geom- 
etry and that only the Ti layers are hydrogenated. The 
agreement between measured and simulated patterns is 
good. Significant discrepancies in the peaks positions 
only occur for the 20x and 40x samples, for which the 
measured peak lies at higher angles with respect to the 
simulated one, indicating that partial release of hydro- 
gen has already occurred from the Ti layers. This is not 
surprising as the 20x and 40x samples have Ti layers as 
thin as 1 and 0.5 nm, respectively. In the simulations 
of Fig. [5] the gaussian distribution fwhm is increased to 
0.6° to account for reduced crystallite size. Significantly, 
while the fwhm value reproduces well the width of the 
fundamental peak for the 40x, 20x, lOx and 5x Mg/TiH2 
multilayers, it is too high for the 2x and lx samples. In 
the latter the coherence length must therefore be larger, 
thanks to the reduced amount of Mg/Ti interfaces. 

The intermediate state, characterized by the loading of 
the Ti layers only, has been explored in detail by measur- 
ing XRR on a lOx sample exposed to ~6 Pa of hydrogen 
at 333 K. Such a low pressure is enough to induce hydro- 
gen absorption in Ti but not in Mg. Figure [5] shows the 
change in electrical resistance occurring upon loading of 
the lOx sample in the in-situ XRR setup at 333 K under 
^6 Pa of hydrogen. The resistance decreases as a confir- 
mation that only the Ti layers are loaded: unlike MgEk, 
which is an insulator, T1H2 is a metal with a higher elec- 
trical conductivity than Tii2£ The temperature transient 
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FIG. 8: (color online) Measured (dots) and simulated (lines) 
XRD patterns for the partially desorbed, "intermediate" 
states of the samples described in the present work. The sim- 
ulations assume samples geometries identical to the original 
multilayers but with fcc-(lll) TiH 2 instead of hcp-(002) Ti. 



highlighted in Fig. [9] is due to the injection of hydrogen 
gas at room temperature, in the sample chamber at 333 
K. Figure [10] shows the XRR pattern measured on the 
lOx sample at 333 K in ~6 Pa of hydrogen. From the 
position of the reflectivity peaks and applying eq. [3] we 
calculate a period expansion of about 5.9%, with respect 
to the as-deposited metallic state (see inset in Fig. [10]) . 
According to the literature values for the molar volumes 
of Ti (V Ti = 10.64 cm 3 /mol) and TiH 2 (V Ti H 2 = 13.3 
cm 3 /mol), the hydrogenation of titanium should expand 
the lattice of the host metal by 25%. Given the 3.6:2 
thickness ratio of Mg and Ti measured by XRR (Table 
ITT]) on the as-deposited sample, a uniaxial vertical expan- 
sion of the period A of 5.9%, due to the hydrogenation 
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FIG. 9: Time evolution of the electrical resistance of sample 
lOx, 10x[Ti(2 nm)Mg(4 nm)]Pd(10 nm), while loading in ~6 
Pa of hydrogen at 333 K in the in-situ XRR setup. 
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FIG. 10: XRR from the partially loaded lOx multilayer at 
Ph 2 = 6 Pa and T = 333 K. The inset shows the measured 
linear dependence of the ratio sin 2 6» n /Ac u kqi versus n . 



of Ti only, implies a vertical expansion of the Ti lay- 
ers of 16.5%. This can only be understood taking into 
account strong plastic deformations and out-of-planc ma- 
terial pile up, due to the hydrogen-induced in-plane co- 
herency stress in the titanium layers *2i This dramatic 
material movement is likely to remove the structural co- 
herence at the Mg/Ti interfaces and to be at the origin 
of the "scissor" effect observed in Mg layers sandwiched 
between Ti thin films 

After measuring the XRR of the intermediate state 
in ~6 Pa of H2 at 333 K, we slowly increased the hy- 
drogen pressure up to 10 3 Pa at constant temperature, 
while measuring the electrical resistance of the film (see 
Fig. [TT]) . The abrupt increase of electrical resistance 
between 100 and 300 Pa is due to hydrogen absorption 
in the Mg layers. After reaching equilibrium in 10 3 Pa 
of H2 we measured the XRR pattern for the fully hy- 
drogenated state, as shown in Fig. [TJ] The fit parame- 
ters used in the simulation in Fig. 112a are given in the 
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FIG. 11: Pressure-resistance-isotherm measured during load- 
ing of the lOx sample, 10x[Ti(2 nm)Mg(4 nm)]Pd(10 nm), at 
333 K in the in-situ XRR setup. 



fourth column of Table [III Upon hydrogenation the bi- 
layer thickness A increases by 23%, from 5.6 nm in the 
as-deposited state to 6.9 nm in the fully hydrogenated 
state. Titanium layers expand by 15%, from 2.0 to 2.3 
nm, in good agreement with the value obtained from the 
period expansion observed upon hydrogenation at 6 Pa 
(16.5%). Magnesium layers expand by ~28%, from 3.6 
to 4.6 nm. In the hydrogenation of bulk Mg the molar 
volume increases by 30%, going from 13.97 cm 3 /mol in 
Mg to 18.2 cm 3 /mol in MgH2: a vertical expansion of 
the Mg layers of 28% is therefore an indication of strong 
plastic deformations. Fig. [13] shows the rocking curves 
measured over the second order reflectivity peak, for the 
as-deposited, intermediate and fully hydrogenated sam- 
ple. The pattern consists of two components, one nar- 
row with a fwhm of 0.023° and one much broader. The 
narrow peak is the specular reflectivity and the width is 
what would be expected from an optically flat sample. 
The broader component comes from off-specular scatter- 
ing which includes scattering from roughness at the inter- 
faces. As can be seen in the figure, the off-specular con- 
tribution increases during loading, indicating an increase 
in interfacial roughnessi 28 i 29 This is consistent with the 
reflectivity simulations which also showed a roughness 
increase upon loading. This roughness increase could 
in principle also be attributed to atomic intcrdiffusion, 
however, since interdiffusion causes a lateral roughness 
with no particular length scale, it contributes zero or a 
constant amount to the off specular scattering. The in- 
crease in the off-specular scattering observed in Fig. [13] 
is therefore due to an increase in roughness in the form 
of thickness variations. The small peaks at either side of 
the specular reflections are likely to be due to off-specular 
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FIG. 12: (color online) XRR pattern of the fully loaded lOx 
sample at ph 2 = 10 3 Pa and T = 333 K. The solid red line is 
a simulation of the multilayer using GenX.— (b) Scattering 
Length Density (SLD) profile corresponding to the simulation 
in Fig. [Hi- 



scattering from the first order reflectivity peaks and not 
to correlations in the roughness. 

Due to the poor crystallinity of the fully loaded sam- 
ples a detailed structural characterization of the hydro- 
genated state cannot be achieved with XRD only. In 
order to measure the crystal phases present in the hydro- 
genated multilayers we performed selected area electron- 
diffraction patterns of a lOx sample covered with only 
1 nm of Pd. Such small amount of Pd allows load- 
ing of the sample but its catalytic activity is rapidly 
suppressed by strong metal-support interaction (SMSI) 
effects^ and it is not enough to promote significant hy- 
drogen desorptionji 3 - The measurements were performed 
both for few seconds and for 10 minutes of electron beam 
exposure. In the former case both tetragonal MgH2 and 
face-centered-cubic TiH2 peaks are observed, while in 
the latter hydrogen desorption from the MgH2 layers 
occurs due to electron irradiation and reflections from 
hexagonal-closed-packed Mg are visible, Fig. [T3] It is 
noteworthy that in the hydrogenated state no signs of 
cubic MgH2 are present: in codeposited Mgj,Tii_ y thin 
films, the hydrogenated state is tetragonal for y > 0.87 
and face-centered-cubic for y < 0.87f^ This fee phase is 
similar to the one proposed for high pressure /S-MgH^ 
and originates from the structurally coherent dispersion 
of Mg-rich and Ti-rich nanosized domains in Mg-Ti thin 
film alloysi^iiS The multilayers studied in the present 
work contain 40 at.% of Ti and one would therefore ex- 
pect a similar cubic structure to occur upon hydrogena- 
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FIG. 13: (color online) Rocking curves measured over the sec- 
ond order reflectivity peak, for the as-deposited, intermediate 
and fully hydrogenated lOx sample. The scans are offset and 
only 50% of the experimental points are shown for clarity. 



tion. Apparently, however, the cubic hydride phase can 
only be stabilized by a very fine 3D dispersion of Mg and 
Ti atoms. A lOx multilayer, consisting of repetitions of 2 
nm of Ti and 4 nm of Mg, is already too "segregated" and 
leads to standard tetragonal MgH2. This result is consis- 
tent with the "scissor" effect observed in Ti-sandwichcd 
Mg film, according to which no clastic interaction exists 
between adjacent Mg and Ti layers^ 
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FIG. 14: Selected area electron diffraction patterns for the 
hydrogenated state of a lOx Mg/Ti multilayer, (a) Short and 
(b) long beam exposure. 



B. Optical spectroscopy 



The loading sequence is further investigated by opti- 
cal spectroscopy measurements. An example is given for 
sample 2x in Fig. [T5j where the measured optical reflec- 
tion of the multilayer, deposited on a quartz substrate 
and covered with 10 nm of Pd, is shown for different 
stages of hydrogenation. Comparison with simulated op- 
tical spectra obtained with SCOUT— shows that the dis- 
appearance of the reflection hump at ~4 eV at the begin- 
ning of the hydrogenation process is due to the formation 
of TiH2 (see white arrow in Fig. \Wk and [T5b). In the 
simulations shown in Fig. [T5]thc thickness increase of the 
individual Ti and Mg layers has been taken according to 
the variations observed by means of XRR, ^15% for Ti 
and 28% for Mg, leading to an excellent agreement with 
the measured spectra. 



C. Hydrogenography 

In order to interpret the Pressure-optical 
Transmission-Isotherms (PTIs) measured by hy- 
drogenography on the Pd-capped Mg/Ti multilayers we 
make use of the following assumptions: 1) Mg layers 
in direct contact with Pd feel an elastic constraint, 
due to the formation of Mg-Pd alloys at the interface, 
which leads to plateau pressures higher than what 
expected from bulk Mg (plf lk K Mg = 12 Pa^)^ 2) in 
the process of Mg-Pd alloy formation typically 6 nm 
of Mg are "lost" and do not contribute to the optical 
change occurring upon hydrogen absorption^ 3) Mg 
films surrounded by Ti layers, on the contrary, do not 
feel significant elastic constraints, thanks to the positive 
enthalpy of mixing of Mg and Ti which leads to poor 
interface adhesion. Furthermore Ti absorbs hydrogen 
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FIG. 15: (color online) Experimental (lines) and simulated 
(dots) reflection spectra measured through the quartz sub- 
strate for the 2x sample, 2x[Ti(10 nm)Mg(20 nm)]Pd(10 nm), 
exposed to slowly increasing pressures of a 6%H2/Ar mixtures 
at room temperature, (a) As-deposited, (b) intermediate and 
(c) hydrogenated state. The dielectric functions of Ti, Mg, 
TiH2, MgEb and Pd used in the simulations are the same 
used in Ref. 14. 



at lower pressures than Mg and the consequent lattice 
expansion removes the partial lattice coherence at the 
Ti/Mg interfaces and leads to quasifree Mg layers 
The PTIs measured at 333 K for the multilayers studied 
in the present work are shown in Fig. 1161 The PTI for 
the lx sample, Ti(20 nm)Mg(40 nm)Pd(10 nm), shows 
a relatively high plateau pressure with respect to bulk 
Mg (Pbuik^g — 12 Pa)r^ due to the elastic clamping of 
the top Pd layer^ The same effect, together with the 
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FIG. 16: (color online) Loading Pressure- Transmission- 
Isotherms measured by hydrogenography at 333 K for the 
Mg/Ti multilayers covered with 10 nm of Pd. 



"scissor" property of Ti^ is responsible for the double 
plateau observed for the 2x sample, Ti(10 nm)Mg(20 
nm)Ti(10 nm)Mg(20 nm)Pd(10 nm): the bottom Mg 
layer (lower plateau) is sandwiched between two Ti 
layers and it is therefore quasifree^ while the top one 
(higher plateau) is in contact with the Pd cover and 
feels its elastic constraint, therefore loading at higher 
hydrogen pressures. Following the same line of reasoning 
one would expect for the 5x sample to find two plateaus: 
a plateau at low pressure, build up by the bottom 4 Mg 
layers sandwiched between Ti and a plateau at higher 
pressure, coming from the top Mg layer in contact with 
Pd. The top plateau however is not visible due to the 
fact that the uppermost Mg layer is only 5 nm thick 
and therefore completely alloyed to Pd. In the lOx 20x 
and 40x the plateau pressure slightly increases with 
respect to the 5x sample. A possible explanation for 
such behavior is that in the lOx, 20x and 40x samples 
the individual Mg layers are very thin (4, 2 and 1 nm, 
respectively) and surface energy might provide a relevant 
contribution to the enthalpy of hydride formation.— 
Another effect that might play a role is the fact that in 
these samples the Ti layers (2, 1 and 0.5 nm) might not 
be thick enough to form perfectly closed films: the Mg 
layers would therefore not be completely shielded from 
the clamping effect of the top Pd cover. For the 40x 
sample a sloping plateau similar to the one measured for 
a co-deposited Mgo.6Tio.4 alloy is observed^ suggesting 
that the microstructure of this sample is closer to a 
mixed alloy than to a well defined multilayer. 

It is noteworthy that for the lOx, 20x and 40x samples 
a small plateau appears at very low hydrogen pressures. 
The width of this plateau is proportional to the number 
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of Mg/Ti interfaces present in the multilayers. We inter- 
pret these plateaus as originating from hydrogen atoms 
populating the interstitial sites located at the Mg/Ti in- 
terfaces: in the metallic films the hydrogen atoms are 
located in interstitial tetrahedral sites, inside the Mg 
layers these are Mg4 sites but, crossing the Mg/Ti in- 
terfaces, there will also be Mg 3 Ti, Mg 2 Ti 2 , MgTi 3 and 
Ti4 sites. Since the formation enthalpy of TiH2 is smaller 
than the one of MgH 2 , the substitution of one (or more) 
Mg atoms with Ti in a tetrahedral site will lower the ab- 
sorption energy and the site will be populated at lower 
hydrogen pressures £ As can be seen in Fig. 111! a simi- 
lar small plateau at the beginning of the hydrogenation 
process is also observed when measuring the film electri- 
cal resistance, while loading the lOx sample used in the 
XRR measurement. A comparison of Fig. QT] and Fig. 
1161 shows that, when loaded in the in-situ XRR setup, 
the lOx sample absorbs hydrogen at a higher pressure 
with respect to the loading in the hydrogenography op- 
tical setup. This is not surprising since, when measuring 
isotherms by hydrogenography, we want to obtain equi- 
librium curves and the pressure is increased very slowly in 
20 hours. In the XRR setup on the other hand the load- 
ing was completed in 30 minutes and it is therefore not 
an equilibrium measurement. The small step observed 
at ~230 Pa in the pressure-resistance-isotherm cannot 
be explained with simple thermodynamic considerations, 
but it could be related to the fact that the measurement is 
done in the kinetic regime, as such a step is not observed 
in the equilibrium isotherm obtained by hydrogenogra- 
phy 



D. Diffusion simulation of hydrogen cycling 

While loading the multilayers in the XRD setup we ex- 
posed them to 1 bar of H2 at room temperature, a pres- 
sure sufficient to hydrogenate both the Ti and the Mg 
layers. Nevertheless we observed an intermediate state 
in which only the Ti layers were loaded. Furthermore, 
upon exposure to air at room temperature hydrogen des- 
orbs from the Mg layers but remains trapped in the Ti 
ones and the original metallic state can only be recovered 
by heating the samples up to 433 K. These loading and 
unloading sequences can be qualitatively explained by 
looking at the difference in thermodynamic properties of 
hydrogen absorption in Mg and Ti. Pasturel et alS^ have 
already shown how the chemical potential of hydrogen in 
different transition metals (TM) can influence the hydro- 
gen sorption kinetics in Mg2Ni/TM/Pd trilayers. Here 
we extend their model to a system with 7 layers, in order 
to account for the multiple repetitions typical of a mul- 
tilayer. The model treats hydrogen dissolved in metals 
in the lattice-gas approximation, without including any 
H-H interaction. The chemical potential of hydrogen in 



a metal (M) is therefore written as: 
^(c M ) = Hh-TSk 

= Hn ~ T (S'conf. + 5 v ibr.) (4) 

Hh and Sh arc the partial molar enthalpy and entropy 
of hydrogen in metals, respectively. The configurational 
contribution to the entropy, S'conf., is proportional to 
In [cm/ (1 — cm)], where Cm is the hydrogen concentration 
in the metal. This term takes care that hydrogen atoms 
obey Fermi-Dirac statistics in the host metallic lattice, 
due to the single occupation of interstitial sites. In most 
metals the vibrational term, SVibr., is small at moderate 
temperatures and can be neglected, with the significant 
exception of Pd (S vihl « 19.2 JK _1 (mol H)" 1 at 298 

Each iteration of the model consists of two steps. First, 
the concentration within each layer is updated according 
to the diffusion equation. Subsequently, we impose the 
equality of chemical potential at the interface between 
two materials by changing the concentrations at the sites 
adjacent to the interface such that the chemical poten- 
tial is made equal on both sites, while keeping the total 
amount of hydrogen atoms unchanged. Note that the 
equality in chemical potentials at the interfaces does not 
imply the equality of hydrogen concentrations in the ad- 
jacent surfaces of the two layers. Instead of using the 
experimental values of the hydrogen pressure and the en- 
thalpies and entropies of hydrogen absorption in Mg, Ti 
and Pd we define the following dimcnsionlcss parameters: 
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where p° = 10 5 Pa, AS^ is the entropy of hydride for- 
mation which, except for Pd (AS£ d = -48.7 JK _1 (mol 
H) -1 ), is taken equal to the entropy of hydrogen gas at 
standard pressure (ASm w — |<$h 2 = ~^5 JK _1 (mol 
H) _1 ) and AHm is the enthalpy of hydride formation: 
AHp d = -20.5 kJ(mol H)" 1 ^ AH Mg = -37.2 kJ(mol 
H)" 1 ^ AH Ti = -65 kJ(mol H)" 1 ^ Substituting these 
values in eq. [5] and taking T = 333 K, we obtain: 
epd = —1.5, CMg = —5.6 and exi = —15.6. 

In Fig. 1171 we simulate the loading and unloading be- 
havior of a multilayer made of 3 Ti/Mg repetitions and 
covered with a Pd layer, in which the Mg layers are twice 
as thick as the Pd and Ti ones: 3x[Ti(z)Mg(2z)]Pd(z). 
In order to highlight only the effect of the chemical po- 
tential, we assume that the hydrogen diffusion coefficient 



11 



is the same for all the materials. The loading pressure is 
taken as 10 5 Pa (en = 0) and the unloading pressure is 
taken negative enough in order to allow desorption from 
all the layers (en = —20). Upon exposure to 10 5 Pa of H2 



Time Time 




Loading Unloading 

FIG. 17: (color online) Simulated time evolution (in loga- 
rithmic scale) of the hydrogen concentration in each layer of 
a 3x[Ti(z)Mg(2z)]Pd(,z) sample, upon hydrogen absorption 
(eH = 0) and desorption (eH = —20). The vertical dashed 
lines represent the intermediate states, observed experimen- 
tally both upon loading and unloading, in which only the Ti 
layers are hydrogenated. The white dash-dot lines are the 
simulated desorption rate from the Ti layers at 433 K. 

the Pd and Ti layers and the uppermost Mg layer load 
very quickly, while the Mg layers "sandwiched" between 
Ti absorb hydrogen more slowly, effectively producing a 
transient intermediate TiH 2 /Mg multilayer. This is due 
to the fact that at the beginning of the loading process, 
when the hydrogen concentration is low everywhere, Ti 
acts as a hydrogen sink due to its lower enthalpy of hy- 
dride formation, effectively sucking hydrogen atoms lying 
at the Ti/Mg interfaces. When the Ti layers are almost 
full the logarithmic entropic term in eq. 2] dominates 
over the energetic term and Mg layers start to absorb 
hydrogen. Upon unloading the effect is opposite: at the 
beginning the hydrogen concentration is very high ev- 
erywhere and hydrogen atoms remain trapped in the Ti 
layers until the adjacent Mg layers are "empty" enough 
to equilibrate the chemical potentials. When desorption 
starts in the Ti layers, however, the hydrogen concen- 
tration in the Mg layers is so low that the total flux of 
hydrogen atoms through Mg (which is proportional to 
the gradient in concentration) is minimal, resulting in a 
very slow hydrogen desorption. In Fig. [T7]wc also show 
that increasing the temperature of the system to 433 K 
enhances the rate of hydrogen desorption from the Ti 



layers (white dash-dot lines). With these simple simula- 
tions we can therefore give a qualitative interpretation of 
the persistence of hydrogen in the Ti layers, both upon 
loading and unloading, as observed experimentally. 



IV. CONCLUSIONS 

We have prepared several Pd-capped Mg/Ti multilay- 
ers with various periodicities by means of magnetron 
sputtering. The deposited samples have lattice param- 
eters in the z direction close to their nominal bulk values 
and low interfacial roughnesses. Partial structural coher- 
ence exists at the Mg/Ti interfaces in the as-deposited 
state but it is reduced upon hydrogen absorption and des- 
orption. The hydrogen loading sequence, as confirmed by 
XRD, XRR and optical spectroscopy, agrees with what 
is expected from thermodynamic considerations on the 
enthalpies of formation of magnesium and titanium hy- 
drides: Mg/Ti -> Mg/TiH 2 -»■ MgH 2 /TiH 2 . Hydrogen 
absorption in both the Ti and Mg layers leads to large 
expansions in the vertical out-of-plane direction, well be- 
yond the elastic regime, indicating that massive material 
pile up, due to plastic deformations and creation of de- 
fects, has to take place. Upon dehydrogenation hydrogen 
is kinctically trapped in the Ti layers. Complete desorp- 
tion only occurs upon exposure to air at 433 K leading 
to the original metallic layered structure, with a shorter 
coherence length. Magnesium hydride in a Mg/Ti multi- 
layer with period A as small as 6 nm, retains its standard 
tetragonal structure, suggesting that the occurrence of a 
cubic hydrogenated phase, as observed in partially seg- 
regated Mg y Tii_ y (y < 0.87) thin filmS) 9 -^ can only be 
stabilized by a very fine dispersion of Mg-rich and Ti- 
rich domains. Pressure-optical Transmission-Isotherms 
measured by hydrogenography, can be interpreted on the 
basis of the clamping effect on thin Mg films due to the 
adjacent layer a 14 i 15 and, possibly, on the surface energy 
differences between the metallic and hydrogenated states 
of ultra thin Mg films. A simple diffusion model allows us 
to reproduce both the loading and unloading sequences 
measured experimentally. 
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